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We present a theoretical analysis of different methods to synthesize entangled states of two super¬ 
conducting resonators. These methods use experimentally demonstrated interactions of resonators 
with artihcial atoms, and offer efficient routes to generate nonclassical states. We analyze the the¬ 
oretical structure of these algorithms and their average performance for arbitrary states and for 
deterministically preparing NOON states. Using a new state synthesis algorithm, we show that 
NOON states can be prepared in a time linear in the desired photon number and without any 
state-selective interactions. 
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I. INTRODUCTION 

In recent years we have witnessed a dramatic evolution in the quantum mechanical experiments performed with su¬ 
perconducting circuits. Initially, the challenge was to fabricate, prepare, and isolate signatures of quantum behavior 
of the coupled motion of Cooper pairs through Josephson junctions and electrodynamic oscillations in supercon¬ 
ducting devices urn- This has now become routine in the field of circuit QED [5], and the frontier is designing 
and manipulating the quantum states of coupled superconducting qubits (quantum bits) and resonators to achieve 
quantum-enhanced information processing laiz]- 

When embarking on this new journey, the quantum mechanical engineer must decide on which degrees of freedom 
does she wish to manipulate: the electronic (qubit) or electromagnetic (resonator)? There are important advantages 
on both sides, but, until recently [8], the coherence (or quality factor) of superconducting resonators (cavities) could 
be significantly greater than the qubit circuits utilizing Josephson junctions. Thus, one should consider what coherent 
operations can be performed with superconducting resonators as opposed to qubit circuits. Such studies include high- 
fidelity measurement EHn, computation [12], and error correction mill], which all attempt to utilize the larger 
state space afforded by the harmonic oscillator states of a resonator to achieve greater efficiency. 

Here we consider how to efficiently manipulate these modes into desired quantum states. In particular, we consider 
theoretical methods to perform “digital” state synthesis of superconducting resonators, where the desired state is a 
superposition of Fock states [T5l[T6]. An alternative “analog” approach uses superpositions of coherent states [T7|. We 
expect that many of the issues encountered in the digital regime will have counterparts in the analog regime, but both 
warrant detailed study. In this paper, we continue the analysis of Fock state manipulations, here for the synthesis of 
entangled states between two resonators. 

The general state synthesis problem concerns how one can prepare, with high fidelity, an arbitrarily chosen quantum 
state. A state synthesis algorithm is a procedure, given a description of the target state, to identify the appropriate 
set of Hamiltonian controls (such as amplitudes and frequencies of control fields) that will prepare the target state 
from a fixed initial state. Note that there are two senses in which the state synthesis problem is solved algorithmically. 
First, a classical algorithm is typically implemented as a computer program to find the set of controls. Second, the 
output of this program is itself a program, namely a sequence of operations to be applied to quantum hardware to 
prepare the desired state. Thus, the state synthesis algorithms presented are a means to program future quantum 
machines. 

In this paper we consider a scenario such as that depicted in Fig. 1, in which a qubit is used to couple two resonant 
cavities A and B, the latter with Fock states \na)0\n}y). We will analyze algorithms that deterministically and exactly 
produce an arbitrarily chosen state of two resonators, of the form 

Na W 

iV^target) ~ |0) qubit ^ E E ^na,ni, \na)<^\nb). (1) 

na=0nb=0 
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FIG. 1: Schematic scenario for entangled state synthesis problems, in which two cavities A and B are coupled by a qubit. The 
qubit can be controlled to effect qubit rotations or swaps between the qubit and each cavity. 


In fact, we will provide a performance analysis of two such algorithms, one a photon subtraction algorithm based 
on previous work miin], and a second photon swapping algorithm new to this work. The results obtained here 
can be used as a benchmark for alternative procedures to prepare such states. These alternatives include numerical 
optimization methods, closed-loop control methods, or other measurement-based methods for state preparation. Fi¬ 
nally, we specifically consider how our algorithm compares with special-purpose NOON-state preparation 1201 EH, and 
demonstrate that our new algorithm can synthesize these and a large class of entangled states without state-selective 
interactions. As state-selective interactions are often weaker than direct qubit-resonator interactions, we expect these 
results will aid future demonstrations of entanglement in superconducting qubit-resonator systems. 

This paper is organized as follows. In Section II, we introduce the general state synthesis problem by studying 
how to prepare a general state of a d-level quantum system (i.e., a qudit). This is followed in Section III by a 
presentation of the Law-Eberly algorithm for a single resonator coupled to a qubit, before addressing in Section IV 
the two entangled-state synthesis algorithms for two resonators coupled by a qubit. Finally, Section V compares 
these algorithms for the preparation of NOON states. We conclude in Section VI by summarizing our work and open 
questions. 


II. QUDIT STATE SYNTHESIS 

Before focusing on state synthesis problems for systems with specific Hamiltonians, it is useful to start with a 
simpler problem. Thus, we begin by considering the synthesis of an arbitrary state of a d-level system known as a 
qudit (a quantum digit) [22]. Such systems are universal for computation [23l [24], and much is known about the 
construction of logic [25H28] and error correction [22l |29] using such systems . For our case, the qudit could be a 
nonlinear oscillator driven directly by control fields with frequencies tuned to distinct transitions, shifted by either 
resonant or dispersive coupling to a qubit [30] . 

Our task is to prepare the quantum state 


d-l 

iV^target) ~ ^ ^ |'^) 

n=0 


( 2 ) 


starting from the initial state |0). For ease of analysis, we write the coefficients of the target state {ip) in terms of d 
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phases and d-dimensional spherical polar coordinates 

Co = cos6>oe*'^° 

Cl = sin6>o cos6>ie*'^^ 


C(i -2 = sin^o sin^i • • • sin ^(^_3 cos 

Cd-i = sin6»osin6»i---sin6»d-3sin6'd_2e*'^‘^-\ (3) 

and the angles have the ranges 0 < Oj < 7r/2 and —tt < (pj < tt. Note that the first phase 0o could be set to zero 
without changing the physical problem. 

The quantum state can be generated by two types of operations: the two-level rotations 


Tln,n+l{d) = exp 


Q 

-i- {\n){n + 1| + |n + l)(n|) 


and the single-level phase shifts 

Znicp) = exp {i(t)\n){n \). 

A solution to the state synthesis problem is the following 

IV^) = Zd-iicpd-i) (Y\^^k{(pk)Zk^i{7T/2)7Zk,k+i{20k) \ |0), 

\k=o / 

where we are using a “time-ordered” product notation, e.g. 

d-l 

l[Uj=Ud-iUd-2---U,. 

i=i 


(4) 

(5) 

(6) 

(7) 


For a qubit (d = 2), this reduces to ^o(0o)>2'i(0i + 7r/2)7^o,i(2do), which can be combined into the product of two 
spin rotations (about the x and ^ axes, respectively). This is the number of rotations required to map an arbitrary 
qubit state’s Bloch vector from the north pole to any point on the Bloch sphere. Similarly, this solution is a minimal 
approach to controlling a qudit, using a fixed set of operations for an arbitrary state of the form Eq. (HI- 

While this solution can be verified by inspection, an alternative approach, the prototype for the state synthesis 
algorithms to be described below, is to find the rotations by reversing the time evolution, that is, to choose a set of 
operations Uj such that 


n = = (8) 

j=d-l 

By simple inversion of Eq. we can use this solution to the inverse evolution equation to find a solution of the state 
synthesis problem given by Eq. The algorithmic approach is to choose each operation to “zero out” an amplitude 
of the target state. Specifically, we index the steps of the algorithm and define the quantum state 




where IpJd) = W and j = d — 1 -> 0. The operator Uj is then chosen so that 

iJl'ipj) = 0. 

Using the rotations specified above, we can set 


where 


aj = avg{{j-ll-ipj+i)) 
Pj = I+arg((j|i/'j+i)) 


= 2 arctan 


O’lV’i+i) 


O' - i|0j+i) 


(9) 

( 10 ) 

( 11 ) 


( 12 ) 
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Before verifying that this produces the same solution as Eq. (|^, let us consider the the first rotation U^_i in Eq. ^ 
(the last rotation of the forward sequence). This is chosen to remove the highest state \d — 1) from the superposition 
in \^d)' Using Eqs. <§ and 0 we have 

{d - IV’) = cos(7d-i/2)cd_i + sin(7d_i/2)cd-2- (13) 

Using the spherical coordinates for q _2 and q_i from Eq. (§ (and cancelling common terms), we thus require 

cos(7^_i/2)sin(9^_2e“*^^-ie*^^-i + i sin(7rf_i/2) cos= 0 (14) 

which is satisfied by 

O^d-l = 

TT 

Pd-1 ~ ^2 ^d-1 

7^_i = 20d-2. (15) 


in complete agreement with Eq. (12). 

The same procedure works for each l/t, the only difference being that the phases of j^i) have already been set 

to zero for j = d — 2, d — 3, • • • (in the previous step), so that in general we find 


Using these angles, we see that 


where 


aj = 

<Pj-i 

Pj = 

TT 

— + dj^d 

7j = 

2d,_i. 

d-l 



\i’) = llUj\o) = Ud-i---Ui\o), 

i=i 


(16) 


(17) 


Uj — (18) 

agrees with Eq. ^ after setting k = j — 1. Note, however, that the choice of the angles is not unique. We could have 
set aj = 0 and pj = 7r/2 + {(j)j — (j)j-i) to achieve the same result. 

Eor an arbitrary target state l^/i), we can characterize the performance of this algorithm in terms of the resources 
needed to construct the state. These resources could be analyzed in terms of the number of controls required, the 
energy associated with each control, and the duration over which the control fields act. Eor simplicity, we consider 
that each two-state rotation can occur with an effective Rabi frequency and each phase shift with TAcj. Then, this 
algorithm produces a set of d — 1 phase shifts (assuming a and (3 can occur in parallel) and d — 1 rotations, such that 
the overall time is 


T 


. d-l . d-1 


i=i 


i=i 


(19) 


The average time required can be found by averaging over the unit circle (for Pj) and the spherical coordinates in Eq. 
^ (for 7 j = 20 j). We find that {/3j) = nj2 and 


Thus, we find that 


{Oj) = 


e dd 

{sin ef~^~^d0 


TT TT 1 

2 “ 4 Vrf - 2 - i ■ 



TT i 


( 20 ) 


( 21 ) 


Thus, this particular sequence takes a time that grows roughly linear in the Hilbert space dimension d, with timescales 
given by l/H and l/Acj. 
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III. LAW-EBERLY ALGORITHM 


Having illustrated the properties of qudit state synthesis, we proceed to a qubit-oscillator system, appropriate for 
superconducting circuits and resonators. This algorithm was first put forward by Law and Eberly in the context of 
cavity-QED [31], and experimentally demonstrated using the internal and vibrational states of a trapped ion [32] . 
The superconducting experiments na HI] demonstrated exquisite control over the combined Hilbert space of the 
qubit-resonator system. 

Here we review this problem, namely how to synthesize an arbitrary state of harmonic oscillator mode (the resonator) 
by using a two-level auxiliary system (qubit). The target state is taken to be 

-/Vmax 

IV’target) = | 0 ) 0 ^ Cn\n), ( 22 ) 

n=0 

in which the resonator has a maximum photon number Amax- The systems are coupled by Jaynes-Cummings-type 
swapping interaction, with a Hamiltonian (in the interaction picture) of the form 

1-Llh = i Acj(t)crV + + g{t) (cr’*‘a -1- aa^) , (23) 

where a = |0)(1| is the lowering operator for the qubit and the control fields (Aco’(t), fl(t), and g{t)) are assumed to 
be under experimental control. These control fields enable unitary operations of the form 

S{0) = exp [—iO (acr^ + aV)] , (24) 


and 


R{0) = exp 



Z(0) = exp 



(25) 


(26) 


The Law-Eberly algorithm will be expressed in terms of these operations. 

The state-synthesis procedure follows a similar pattern as the qudit case presented above. We first set 

|V’i) = C/j|t/'i+i) (27) 

where 

uj = R\^j)ZH(3j)S\9^)ZHaj). (28) 

and IV^AT+i) = IV^target)- Here and 0 are chosen so that at each step 

(0,j#i) = (l,i|^j)=0, (29) 


These angles are then found for each j = N ^ 1, after which = |0,0). The inverse sequence specifies how to 
prepare the target state using only qubit rotations, phase shifts, or qubit-resonator swaps. 

To see how this can be accomplished, it is convenient to break Eq. into two steps by defining 

l^j+ 1 / 2 ) = S\d^)Z\a^)\^^+,), (30) 

and 

IV’i) = RHlj)zHl3jMj+i/2)- (31) 


Eor convenience, we also define 


i’q.kU) = («> where g = 0 or 1. 


(32) 


The first step solves ipojij + 1/2) = 0. Using Eq. (30), this reduces to 


cos (v^6»j) ?/)oj(j + 1) + ze sin (v^6»j) -ipij-iij + 1) = 0 


(33) 
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or 


V V’i,i-i(j + l) 


This is solved by 


aj = arg 


Oj = arctan 

vi 


(i,j - i|V>j+i) \ 

(ojli’j+i) 


{i,j - i|t/>j+i) 

The second step solves = 0. Using Eq. ( [M] ), this reduces to 

cos ipij-iU + 1/2) + sin i’oj-iU + 1/2) 


or 


tan ( ^ 


(I) 


This, in turn, has the solution 


j3j = arg 


7 j = 2 arctan 


+ 1/2) 
V’o,j-i(j + 1/2)’ 

z(l,j-l|^j+l/2) \ 

(0,i - l|i/’j+i/2) / 
(l,i-l|^j+i/2) 


(0,i - l|i/’j+i/2) 


(34) 


(35) 


(36) 


(37) 


(38) 


By solving these equations for aj^/3j^jj, and Oj for j = N ^ 0, keeping track of l^pj) at each step, the amplitude 
is forced down to smaller and smaller photon numbers, so that l^po) = |0,0). The form of the sequence was chosen 
specifically to not send amplitude to higher photon numbers. That is, the algorithm actually produces operators Uj 
and states \ 2 pj) that satisfy the condition 


(0, kl'ipj) = (1, kl'ipj) = 0 for k > j. 


(39) 


This condition is the most challenging to generalize to more resonators. 

The average values of and Oj can be used to find the average time required, assuming constant controls 

TAcj, O and g: 


-I -^max -j -^max A/max 

^ ^ ^ S + l^il) + o S u + - ^ 


i=i 




(40) 


i=i 


i=i 


We illustrate the average angles obtained with Eqs. (35) and (38) in Fig. Here we have generated one hundred 
random target states for each value of Amax and averaged the total of the angles used in the Law-Eberly algorithm. 
Also shown are the approximations 


j 

S TT ^A^max 2 ^ ’ 


j 

s 2.72Af„,ax - 1.66, 



2.65v/y//; - 1.78, 

(41) 


3 


obtained by fitting the numerical data. 

The linear increase of the phases and qubit rotations are expected, as each step requires such a rotation, while 
the square-root dependence of is due to the y^i-coupling between the qubit and the n-photon state of the 

resonator. When these averaged angles are substituted into Eq. (40), we see that, just as the qudit case, an arbitrary 
state of the form Eq. (22) can be synthesized in a time proportional to the effective Hilbert-space dimension. 
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FIG. 2: Averaged total angles for the state synthesis sequence using the Law-Eberly algorithm. 


IV. TWO-RESONATOR ALGORITHMS 

The state synthesis problem can be extended to any number of resonators, but explicit algorithms are a challenge 
to specify. Early work utilized special interactions [33H37] to enable the transfer of excitations between resonators 
and multi-level atoms. These interactions, while natural to trapped-ion systems, are not directly applicable to the 
Hamiltonians considered here. There have been a number of recent studies of interesting interactions that can be 
generated between superconducting or nanomechanical resonators [38ti4^ . As we are interested in Fock-state control, 
we consider the simplest system of a single qubit coupling two resonators, with each interaction of the Jaynes- 
Cummings form. We further consider algorithms that accomplish the synthesis of an arbitrary two-resonator state 

Na A^b 

iV^target) = |0) G) E E Cna,n{, \na) 0 \nb). (42) 

n„=0 nf,=0 

In general, such a state will be entangled, thus we call this the entangled-state synthesis problem. 

The interactions used in these algorithms are all based on the underlying Hamiltonian: 

H/h = -b + ga{t) {cr^a + aa^) + + ab^) . (43) 

For our purposes, we again assume that the control fields (Aco’(t), ga{t)^ and gb{t)) can be turned on and off at 
will, and thus we restrict our attention to swap operators 

A{0) = exp [—iO (acr^ -f aV)] , 

B{0) = exp [—iO (ba^ + ^ 

the single-qubit phase rotations 

Z{(p) = exp ) (45) 

and the number-state-selective qubit rotations 

Rn^,nt(.^) = exp f -■* 2 '^® 0 |na,nb)(na,nb| 

The last operation utilizes the Stark-shift of each resonator on the qubit, and can (in principle) be extended to 
many resonators. The actual operation may or may not be completely selective on an individual Fock state 
but the basic conditions required can be satisfied provided there is some selectivity. Such selective operations were 

first observed in circuit QED as number splitting [43] and later used for photon measurement [9] and theoretically 

proposed for state synthesis m- The detailed physics of such interactions were further analyzed for both qudit 
operations [12] and for state synthesis m- We will assume complete selectively here, but will discuss how entangled 
states can be synthesized with reduced selectivity in the next section. 

We will be using a Fock-state diagram, such as Fig. in which a state \q^na^nb) with Ua excitations in mode A, 
rib excitations in mode H, and qubit state q is indicated by the node at location {ria^Tib) and internal level q. Each 


(46) 


(44) 













































of the operations described above corresponds to a transition between sets of states in this diagram, and the state 
sythesis sequence can be interpreted using paths in this diagram. Two algorithms, to be described below, can be 
visualized using these diagrams. The first algorithm, which we call the photon substraction algorithm, uses vertical 
and horizontal paths from top-to-bottom and left-to-right in the Fock-state diagram. Thee second algorithm, which 
we call the photon swapping algorithm, uses diagonal paths from the upper-left to the lower-right. These are the 
two natural choices for how to navigate the Fock-state diagram in order to program the quantum system into any 
desired state. In this section, we will analyze each algorithm in detail, and compare their average performance when 
preparing an arbitrary two-resonator state of the form Eq. (42). 


A. Algorithm 1: Photon Subtraction 

The first algorithm for superconducting resonators m used a strategy similar to the trapped-ion proposal by Kneer 
and Law [36], and involves repeated subtraction of photons from one of the resonators. Using the Fock state diagram 
presented in Fig. state amplitudes are cleared column-by-column, row-by-row, until all of the remaining photons 
are in one mode only. The final steps remove these photons by the Law-Eberly protocol described above. 

The essential steps can be written as 


u = Ua (47) 

where 

Na 

Ua = ll Ziaj)A{ej)Zil3^)Rij^) (48) 


and 


Nb 


^{^3k)B{0jk)Z{Pjk)Rna=k{ljk) 


(49) 


/c=0 


Read in reverse, the elements of Ua and are all of the form of Eq. (28), with phases and angles calculated 


using the same method. In more detail, Ulj is a product of operations that subtract a photon from state |0,/c,j) 
(transferring its amplitude to |1, /c — 1, j) and then to |0, /c — 1, j)), first for k = Na ^ I (column-by-column), which is 
then repeated for j = ^ 1 (row-by-row). After all of the amplitudes have been transferred to the states |0, k,0), 

Ua removes these much as the original Law-Eberly algorithm. A graphical representation of this sequence is shown 

in Fig. 1^ 

To prevent amplitudes from returning to previously cleared states, it was proposed [18] to make the qubit rotations 
in number-state selective. This could be achieved by choosing a rotation for state (n^ = k^niy = j — 1) only, 
but the main requirement is that previously removed states with rib = j and ria < k are unaffected. Note also that, 
for certain types of number-state-selective interactions, the column ordering may need to be reversed, as discussed in 

m- 

The actual steps involved in this algorithm are nearly identical to those in the Law-Eberly algorithm. The main 
challenge is to keep track of the various quantum states and the ordering of the operations. Eor completeness, we 
include an explicit treatment here, first breaking up the quantum evolution into two stages (for the B and A swaps, 
respectively). Eor the first stage, we define 


l^j,fe+i/2) = B^ejk)Z\ajk)\ipj ,/c + l) 

1^4^j,k) = l^i,^+i/2) 


(50) 


where k is the “fast” index (ranging from Na 0) and j is the “slow” index (ranging from N^ ^ 1). These states 
have the boundary conditions |'0j^7Va+i) = |'0i+i,o) and \'ipNb,Na-\-i) — IV^target)- Eollowing a procedure similar to the 
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FIG. 3: Illustration of the state synthesis sequence using the photon subtraction algorithm. 


previous section, we find 


^jk 


Ojk 

Pjk 

Ijk 


arg 


arctan 

VJ 


i{0,k,j\^pj^k+i) 


arg 


2 arctan 


{l,k,j - l\'ipj,k+i) 
i{l,k,j - l|t/>j,fc+i/2) \ 

{0,k,j - l\'ipj,k+i/2) J 

{l,k,j - l\i’j,k+i/2) 


{0,k,j - l\'ipj,k+i/2) 


These equations can be solved for k = Na ^ 0, j = ^ 1, until we reach the second stage. 

For stage two, we define 


(51) 


lV’i+ 1 / 2 ) = A'<{ej)Z\aj)\ipj+i), 

It/.,-) = (52) 

with j ranging from Na 1 and |'07v^+i) = (the final state from stage 1). The remaining parameters are then 

found by 


= 




1 


On = arctan 


Vj 


{0,j,0\'ipj+i) 


{l,j - l,0\'ipj+i) 


Pi 

7i 


(OJ^WTSW)) 


= 2 arctan 


(i,i - 1 , 01 -^,-+1/2) 


(o,i -1,01-^,-+1/2) 


(53) 


The total number of operations amounts to Na-\-Ni,-\-NaNi, swaps, Na+Na+NaNi, rotations, and 2{Na-\-Ni,-\-NaNi,) 
phase shifts. Assuming we can turn the various Hamiltonians on and off with rates TAcj, and (for the phase. 
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FIG. 4: Averaged total angles for the state synthesis sequence using the photon subtraction algorithm. 


swap, and rotation operators, respectively), the total time for this sequence is 


T = 




Auu 


jk 


5”)-K 




jk 


(54) 


jk 


The averaged total angles are shown in Fig. These were again formed by generating one hundred random target 
states of the form Eq. (42) with Na = Nij = A^max and summing and averaging the angles produced by Eqs. (51) and 
(53). Also shown are the approximations 


n 

X(^") 

n 

X(^") 


2.87VL, + 4.67V^ax-2.7, 
2.97Vmax - 0.7, 


(55) 


where the sums are over all of the indices (j, k) and j of stages one and two, respectively. The quadratic growth of the 
phase shifts and rotations match the total number of operations, while the linear growth of the swap angles is reduced 
by a square-root, similar to the Law-Eberly results above. We observe that the total time required to produce an 
arbitrary target state again scales with effective Hilbert-space dimension. 


B. Algorithm 2: Photon Swapping 

An important fact regarding state synthesis is that the solution need not be unique. There are an infinite number 
of solutions, and finding the optimal solution, under constraints on time, energy, or complexity, is a hard problem. 
Eor the algorithm just presented, we observe that each qubit rotation can add or remove one quantum of energy, 
which occurs A^^ax tiffies. The desired state, however, has a maximum energy of 2Amax (for the resonator state 
|A^max) G) |A^max))- Thus, wc might expect an optimal solution would have a smaller number of qubit rotations. We 
will provide such a solution in this section. Eor convenience, we will let the target state be a superposition of all states 
with quantum number < 2Arnax- 

l^) = |0)® X Cn^,nMa) ®\ni,). (56) 

The photon subtraction algorithm attempts to remove energy for each and every possible state. However, one can 
just as easily move about the Eock state diagram by swapping photons between the resonators, either directly or 
through the qubit [7]. Using the latter, we can swap photons along diagonal paths with a fixed number of quanta 
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FIG. 5: Illustration of the state synthesis sequence using the photon swapping algorithm. 


{ua + + g). By first swapping all of the photons from resonator b to a, we can then remove one photon at a time, 

using only 2A^rnax qnbit rotations. Specifically, our new algorithm is 

2Vniax 

U= U Ue, (57) 

£=1 


where 


1 

m=i 


(58) 


The interpretation of each operation is analogous to Algorithm 1, however, the sequence of operations is significantly 
different. A graphical illustration of this photon swapping algorithm is presented in Fig. 

In this algorithm, as we move along the diagonal path with + g the operator {r]rn,£-m) implements the 

transition 10, ^ — m, m) ^ |1, ^ — m, m —1), while for m > 1 implements the transition |1, ^ — m, m —1) ^ 

|0, ^ — m + 1, m — 1). This sequence has the effect of repeatedly swapping quanta from mode B to mode A, until we 
reach |1,^ — 1,0). At this point, 6>o/-i is chosen to complete the swapping transition |0,-^,0) ^ |1,^ — 1,0), which is 
finally rotated to |0,^ — 1, 0). This last step need only be selective on = ^ — 1, and is so indicated in Uj. 

For completeness, we present the detailed steps of the algorithm. We again break each step in two by defining 

\'^£,m-\-l/2) ~ (^m/—m) |'0-^,m+l) 

and the various angles are calculated by the following equations: 


^m,£—m — arg 

1 

Vm,£—m )F= 

Vn 

Pm—i,£—m — arg 
^m—l,£—m 


i{0,£- m,m\'ipe,m+i) 


arctan 


(0,^ - m,m\'tpe^rn+i) 


{l,£-m,m-l\'ipe,m+i) 


i{l,e - m,m - l\ipe^rn+i/2) 

{0,i-m+l,m- l\ipe,m+i/ 2 } 


\/£ — m + I 


arctan 




{0J-m+l,m- l\ipe,m+i/ 2 ) 


(60) 
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FIG. 6: Averaged total angles for the state synthesis sequence using the photon swapping algorithm. 


These can be solved from m = £ until m = 1, for which we must modify our equations by 


^o/-i 

^o/-i 


arg 


^ arctan 

VI 


{l,£-l,0\^e,i+i/2) \ 
1+1/2) ) 

(0, £, 0|V’f, 1 + 1 / 2 ) 


{iJ-iAV,i+i/2) 


This still leaves a final phase and amplitude rotation, the latter selective on — 1, with parameters 


(61) 


li 


arg 


\ ( 0 ,^-i, 0 |V’«,i) / 


2 arctan 


(i,^-i,o|V>^,i) \ 
(o,^-i,o|V’^,i) ) 


(62) 


This sequence is then repeated for the next diagonal with Ua + rn, + q = ^ — 1, starting with m = £ — 1 and 
Ve-iV = ViV, and again for £ = -+ 1. 

The expectation values for the sum of the angles, when averaged over many target states, are shown in Fig. 
along with the approximate forms 


X](l«ra| + l^ral + \4>n\) 
n 

n 

^ Vn) 

n 


6 . 4 iV^ax, 
bA^max - 3, 

1.657V2 +4.27V^ax-4.5. 


(63) 


Here we see that the total rotation angles is now linear with the maximum photon number, at the cost of an 

increased number of swaps (and phase rotations). However, this can represent a significant advantage, as the qubit 
rotations have (so far) been required to be number-state-selective, and thus limited in Rabi amplitude ft [19]. By 
reducing the number of such rotations, the total time can be reduced. The overall scaling of the time, however, is 
again proportional to the effective Hilbert-space dimension. 


V. NOON-STATE SYNTHESIS 

A natural target state to characterize entangled-state synthesis algorithms is the so-called NOON state 

l^target) = |0) ® ^ d V 0) + |0, N)) , (64) 

an entangled superposition of resonator states in which N photons are in mode A or mode B. This state can be 
considered a generalization of the Bell and GHZ states, and has potential applications in quantum metrology [44] . 
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(a) (b) 




FIG. 7: NOON state synthesis (with iV = 3) using (a) the photon subtraction algorithm, and (b) the photon swapping 
algorithm. Both prepare the NOON state with a linear number of steps, but the latter uses no state-selective interactions. 


An initial experiment to generate a ‘high’ NOON state (with N = 3) was performed using a particular preparation 
method [20l[2T]. That method uses a pair of three-level systems to couple the resonators, and a state-selective swap 
which limits the coupling rate due to the anharmonicity of these system [12] (refinements of this method [45j [46| 
also use state-selective swaps). Comparison of that approach with the original state-synthesis algorithm [18] (the 
photon subtraction method presented above) showed that both were experimentally comparable as far as decoherence 
is concerned [19] . Here we consider the new state-synthesis algorithm presented here and show that no state-selective 
interactions are required, allowing for faster operations and simplified experimental design. Furthermore, using this 
photon swapping algorithm any state of the form 


N 

iV^target) = |0) G) ^ ^ | A" Tl^Tl) 

n=0 


(65) 


can be synthesized without state-selective interactions. 

This improved performance is due to the paths through the Fock-space diagram taken by the new algorithm. By 
following the photon swapping method, starting from a superposition of a given diagonal = A, the first time 

through one can move all of the population down to |0, A — 1,0). Thus, one need only use a Law-Eberly sequence 
along the path = 0 to remove the photons from the system. The result is that any “diagonal” state of the form 
Eq. (65) can be synthesized by one sequence of photon swaps followed by a Law-Eberly sequence with no state- 
selective interactions. The specific set of parameters for NOON state synthesis with A = 3 are shown for the photon 
subtraction and swapping algorithms are shown in Tables [Tj and [Til respectively, and graphically represented in Eig. 

Ill 


We now compare these two approaches for a general NOON state. Based on previous analysis [19], we find that 
the photon subtraction method requires A A-swaps, A H-swaps, and 2A rotations. No phase shifts are required, and 
the parameters scale as 



subtraction 


subtraction 


TT 




6 a/ Amax 3.3. 


( 66 ) 


As in the general photon subtraction algorithm, each of the rotations must be number-state-selective. 

The photon swapping algorithm requires 2A — 1 A-swaps, A H-swaps, and A rotations. There are also a few phase 
shifts required, but they do no scale appreciably with A. By looking at the numerical performance for the NOON 
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TABLE I: NOON State Synthesis by the Photon Subtraction Algorithm 


Step Parameters Quantum State 



71 = 7Tl2,na = 0,726 = 0 

| 0 , 0 , 0 )-i 

1 , 0 , 0 ) 

^1 

II 

to 

| 0 , 0 , 0 )- 

0 , 1 , 0 ) 

R 2 

72 , = 1,726 = 0 = TT 

| 0 , 0 , 0 ) + z 

1 , 1 , 0 ) 

A 2 

02 = 7r/2V2 

| 0 , 0 , 0 ) + 

0 , 2 , 0 ) 

Rs 

73 = TT, 72a = 2, 726 = 0 

|0,0,0>-7 

1,2,0) 

A 3 

03 = 7r/2V3 

|0,0,0)- 

0,3,0) 

Ra 

74 = 7r,72a = 0,726 = 0 

-i|l,0,0)- 

0,3,0) 

Bi 

04 = n/2 

-|0,0,1)- 

0,3,0) 

Rs 

75 = TT, no = 0 , 716 = 1 

7|1,0,1)- 

0,3,0) 

B 2 

05 = n/2V2 

|0,0,2)- 

0,3,0) 

Rs 

76 = TT, no = 0 , n 6 = 2 

-i|l,0,2)- 

0,3,0) 

Bs 

03 = 7r/2V3 

-|0,0,3)- 

0,3,0) 


TABLE 11: NOON State Synthesis by the Photon Swapping Algorithm 
Step Parameters Quantum State 

Ri 71 =^ -^|1,0,0) 

Ai ^i=7r/2 -|0,1,0) 

R 2 72 = tt +^| 1 , 1 , 0 ) 

A 2 02 = 7r/2^/2 +|0,2,0) 

Rs 73 = TT -i| 1,2,0) 

As 03 = 0.2153 -0.3643|0,3,0>-10.9313|1,2,0> 

Bi r?i = 2.1999 -0.3643|0,3,0) + 10.548|1, 2, 0) - 0.753|0,2,1) 

A 4 04 = 1.3589 +0.6454|0,3,0) - i0.1283|l, 2, 0) + 0.2589|0, 2,1) + i0.7071|l, 1,1) 

B 2 r )2 = 7r/2v^ +0.6454|0,3,0) - i0.2889|l, 2, 0) + 0.7071|0,1, 2) 

As 05 = 7r/2 -0.7071|0,3,0)-10.7071|1,0,2> 

Bs r )3 = -0.7071|0,3,0) - 0.7071|0,0, 3) 


state (not shown), we find 


« -kN 

swapping 

« 9.9x7]^-9.2 (67) 

swapping 

As described above, these rotations need not be number-state-selective. 

Comparing these two algorithms for the NOON state, we see that we have traded rotations for swaps, with the 
photon swapping algorithm achieving the optimal number of rotations. Even better, the photon swapping method does 
not require those rotations to be state-selective. Thus, the method presented here has advantages both theoretically 
and experimentally, with the potential for fast performance and optimal scaling. 



VI. CONCLUSION 

In this paper, we have studied state synthesis algorithms for superconducting resonators. By reviewing the qudit 
and Law-Eberly schemes, we have shown how solving for the inverse evolution allows one to determine the operations 
needed to synthesize an arbitrary state. We have further shown how these step-by-step procedures have a complexity 
that typically grows linearly with the effective Hilbert-space dimension. These schemes have been extended to two 
different state synthesis algorithms for a qubit coupled to two resonators. The first type uses photon subtraction 
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to ensure that the inverse evolution leads to the ground state, whereas the second uses photon swapping before any 
photons are removed from the system. When taken in reverse, these algorithms allow one to synthesize an arbitrary 
entangled state of two resonators. Finally, when applied to typical superconducting circuit experiments, we expect that 
the photon swapping method will have improved performance due to a reduced number of state-selective interactions. 

While we have found an improved algorithm, we cannot claim to have found an optimal algorithm. Indeed, we have 
reason to believe that numerical optimizations using the same basic Hamiltonians can lead to improved methods for 
state synthesis. However, we also have reason to believe that the two algorithms compared here are the most natural 
analytical approaches to state synthesis. At the same time, the differences between the two algorithms suggest that 
different types of optimizations may be possible. The photon swapping algorithm minimizes the number of A-swaps 
performed on the system, but at the cost of a quadratic number of H-swaps and state-selective qubit rotations. By 
constrast, the photon swapping algorithm minimizes the number of rotations, at the cost of an increased number 
of A-swaps and slightly increased overall complexity. Nevertheless, for states such as the NOON state, the photon 
swapping method appears to have overall better performance, in that no state-selective rotations are needed at all. 

Finally, the linear scaling of the NOON state sequences are nearly ideal, in that the energy of the final state and 
the number of qubit rotations used (to put energy into the system) are both linear in the state number A^max [E] • We 
further observe that one can achieve a reduction in time complexity by a factor of two by driving multiple transitions 
simultaneously m, but the linear scaling remains. However, recent work has found, using numerical optimization, 
that sublinear scaling is possible for Fock state preparation by starting from a large-amplitude coherent state and 
optimized displacements of the resonator [47]; extending such a scheme to NOON state synthesis is an interesting 
question. 

In conclusion, we have improved the theoretical understanding and performance of entangled-state synthesis algo¬ 
rithms for superconducting resonators. We hope that the results presented here, on a fundamental quantum control 
problem, may provide useful benchmarks for future explorations of control of superconducting or other resonator-based 
systems. 
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